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Any symmetry reduces a second-order differential equation to a first integral: variational sym- 
metries of the action (exemplified by central field dynamics) lead to conservation laws, but symme- 
tries of only the equations of motion (exemplified by scale-invariant hydrostatics) yield first-order 
non-conservation laws between invariants. We obtain these non-conservation laws by extending 
Noether's Theorem to non- variational symmetries and present an innovative variational formulation 
of spherical adiabatic hydrostatics. For the scale-invariant case, this novel synthesis of group theory, 
hydrostatics, and astrophysics allows us to recover all the known properties of polytropes and define 
a core radius, inside which polytropes of index n share a common core mass density structure, and 
outside of which their envelopes differ. The Emden solutions (regular solutions of the Lane-Emden 
equation) are obtained, along with useful approximations. An appendix discusses the n = 3 poly- 
trope in order to emphasize how the same mechanical structure allows different thermal structures 
| ^ ' ■ in relativistic degenerate white dwarfs and zero age main sequence stars. 

\ PACS numbers: 45.20.Jj, 45.50.-j, 95.30.Lz, 97.10. Cv 
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I. SYMMETRIES OF DIFFERENTIAL EQUATIONS AND REDUCTION OF ORDER 

c3 

Noether's Theorem relates every variational symmetry, a symmetry of an action or similar integral, to a conservation 
law, a first integral of the equations of motion By an extension of Noether's Theorem, non-variational symmetries, 
— symmetries of the equations of motion which are not in general variational symmetries — also lead to first integrals, 
which are not conservation laws of the usual divergence form, as discussed in a previous article Q • There it was shown 
that a Lagrangian L(t, qi, q\) and action S = J C(t, q%, 4i)dt, with degrees of freedom q^, can be transformed under an 
infinitesimal point transformation 5(t,qi),Sqj(t,qi): 

_\dG r d(St) 
dt Hl dt ' " 

in terms of the total derivative of the Noether charge, G :— C ■ St + pi ■ [6qi — q\5t), and the variational derivative 
T>i := dC/dqi — d(dC/ 'dq\) / 'dt. For transformations that leave initial and final states unchanged, the variation in 
action is 



(N 
(N 



SC = CSt + (dC/d qi )S qi + (dC/dq,) - q~] = - £ • ^ + 2? 4 • (<% - q t 6t) . ( 1 1 



X 



/ f r / dW BC \ i 

5S if = G(f)-G(i)+ dt[s qi -V i + St-^ + ^\ , (2) 



if the term in d(St)/dt is integrated by parts. If the system evolution obeys an action principle, that this variation 
vanish for independent variations Sqi,5t that vanish at initial and final times, the system obeys the Euler-Lagrange 
equations Oj = and dH/dt = —dC/dt, the rate of change of the Hamiltonian in non-conservative systems. On-shell, 
where Dj = 0, 

SS if = j ' SC dt = G(f) - G(i) (3) 
dC 

— =6C:=SC + C- (dSt/dt). (4) 
dt 

This is Noether's equation, giving the evolution of a symmetry generator or Noether charge, in terms of the Lagrangian 
transformation that it generates. It expresses the Euler-Lagrange equations of motion as the divergence of the Noether 
charge. This divergence vanishes for a variational symmetry, but not for any other symmetry transformation. 
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Noether's equation (4) could have been derived directly from the definition of the Noether charge. But using the 
action principle make manifest the connection between Noether's equation and the Euler-Lagrange equations. We 
use the action principle and this connection to reformulate the theory of hydrostatic barotropic spheres, which is 
integrable if they are scale symmetric, even where this scale symmetry is not a symmetry of the action (Section II). 
The first integrals implied by any symmetry of the equations of motion, while generally not vanishing-divergence 
conservation laws, are still useful dynamical or structural first-order relationships. 

Because it neglects all other structural features, scaling symmetry is the most general simplification that one can 
make for any dynamical system. For the radial scaling transformations we consider, 5r = r, the Lagrangian scales as 
some scalar density 8C — —2ui£; and the action scales as SS = (1 — 2Q)S. The Noether charge generating the scale 
transformation evolves according to a non- conservation law dG/dt = (1 — 2Q)C, a first-order equation encapsulating all 
of the consequences of scaling symmetry Q . From this first-order equation follow directly all the properties of index- n 
polytropes, as established in classical works [!, H|, modern textbooks [1, Q, and the recent, excellent treatments of 
Horedt and Liu 

Our secondary purpose is to present an original variational formulation of spherical hydrostatics and to extend 
Noether's Theorem to non- variational scaling symmetry, which yields a scaling non- conservation law (Section II). For 
spherical hydrostasis, we define a core radius, inside which all stars exhibit a common mass density structure. Outside 
this core, polytropes of different index n show different density structures as the outer boundary is felt (Section III). 
Section IV completes the integration of the Lane-Emden equation by quadratures and obtains useful approximations 
to the Emden function n (£). 

An appendix reviews the thermodynamic properties of the physically important polytropes of index n = 3 @, d, @| . 
What is original here is the explanation of the the differences between relativistic degenerate white dwarf stars and 
ideal gas stars on the zero-age main sequence (ZAMS), following from their different entropy structures. Our original 
approximations to 9a(£) should prove useful in such stars. 



II. SCALING SYMMETRY AND INTEGRABILITY OF HYDROSTATIC SPHERES 



A. Variational Principle for Hydrostatic Spheres 



A non-rotating gaseous sphere in hydrostatic equilibrium obeys the equations of hydrostatic equilibrium and mass 
continuity 

— dP/pdr = Gm/r 2 , dmjdr = ^r 2 p , (5) 

where the pressure, mass density, and included mass P{r), p(r), m(r) depend on radius r. For dependent variables, 
we use the gravitational potential V(r) = J Gm/r 2 dr and the thermodynamic potential (specific enthalpy, ejection 

energy) H(r) — dP/p, so that ([5]) and its integrated form become 

- dH/dr = dV/dr, V{r) + H(r) = — ^ , (6) 

expressing the conservation of the specific energy as the sum of gravitational and internal energies, in a star of mass M 
and radius R. The two first-order equations ([5]) are equivalent to a second-order equation of hydrostatic equilibrium, 
Poisson's Law in terms of the enthalpy H(r): 

?|(- , f)+<"<W>-° ■ < 7 » 

We assume a chemically homogeneous spherical structure, and thermal equilibrium in each mass shell, so that 
p(r), P(r), H(r) are even functions of the radius r. At the origin, spherical symmetry requires dP/dr = and mass 
continuity requires, to order r 2 , 

p(r)=p c (l-Ar 2 ), m (r) = ^f-(l-^Ar 2 ) = ^f.p 2 J 5 p(r) 3 / 5 . (8) 

The average mass density inside radius r is p(r) := m(r) / '(47rr 3 / '3) = p 2 J^ p(r) 3 / 5 . 

In a previous paper Q , we showed that hydrostatic equilibrium (J7J follows from the variational principle SW = 
minimizing the Gibbs free energy, the integral of the Lagrangian C : 

W := / drC{r,H,H') , (9) 
Jo 
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£(r,H,H')=4:7rr 2 [-H /2 /8irG + P{p)}dr , ' := d/dr , (10) 

W is the sum of the gravitational and internal specific energies per radial shell dr. The canonical momentum and 
Hamiltonian, 

to := dC/dH' = -r 2 H'/G , H(r, H, to) = -Gm 2 /2r 2 - 4irr 2 P(H) , (11) 

are the included mass and energy per mass shell. The canonical equations are 

dH/dm = H' = -Gm/r 2 , dU/dH = -m = -Anr 2 p . (12) 

Spherical geometry makes the system nonautonomous, so that dH/dr — —dC/dr — —2C/r vanishes only asymptoti- 
cally, as the mass shells approach planarity. 

The equations of hydrostatic equilibrium (5) can be rewritten 

d\ogu/d\ogr — 3 — u(r) — n(r)v(r) , d\ogv/d\ogr = u — 1 + v(r) — d\og [1 + n{r)]/d\ogr , (13) 

in terms of the logarithmic derivatives 

u{r) := d\ogm/d\ogr, v(r) := —d log (P/p)/d log r, w(r) :— n(r)v(r) — — dlogp/dlogr , (14) 

and an index n(r) 

n{r):=d\ogp/dlag(P/p) , 1 + -^r := dlogP/dlogp , (15) 

n(r) 

which depends on the local thermal structure. The mass density invariant w makes explicit the universal mass density 
structure of all stellar cores, which is not apparent in the conventional pressure invariant v. 

B. Scaling Symmetry and Reduction to First-Order Equation Between Scale Invariants 

Following the our results [2J, a hydrostatic structure is completely integrable, if the structural equations (5) are 
invariant under the infinitesimal scaling transformation 

Sr = r, 5p=—n£j n p, SH = — Q n H, 6H = — (1 + 6j n )H , where uj n :— 2/(n — 1) , (16) 

generated by the Noether charge, for constant n, 

G n :=-U-r-m- (Q n H) = r 2 [(^ + 4nP(H)) ■ r + ■ (17) 

The Lagrangian (10) then transforms as a scalar density of weight — 2£j n 

5C = -2u n C , SS 12 = (1 - 2u n ) ■ S 12 , (18) 

so that only for the n — 5 polytrope (Q n — 1/2) is the action invariant and scaling a symmetry of the action. 

Both structural equations (fTB")) are autonomous, if and only if n is constant, so that, P(r) = K p(r) 1+1 / n , with the 
same constant K (related to the entropy) at each radius. When this is so 

du/d log r — u(3 — u — w n ) , dw n /d\ogr = w n (u — 1 + w n /n) (19) 

d\ogw n d\ogu dlogm 

— dlogr = . (20) 



u—l + Wn/n (3 — u — w n ) 

In this section, we consider only the first equality in (20) 

dw n _ w n (u-l + w n /n) 
du u{3 — u — w n ) 



(21) 



between scale invariants, which encapsulates all the effects of scale invariance. We consider only simple polytropes 



with finite central density p c , so that the regularity condition ([5]) requires that all w n (u) be tangent to | (3 — u) 
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at the origin. Such Emden polytropes are the regular solutions w n (u) of the first-order equation (19), for which 
w n {u) — > | (3 — u) for ii — y 3. 

In terms of the dimensional constant, dimensional radius, and the central enthalpy and pressure 

c? := ij ^r K Pl /n - X > r:=a£, H c := (n + l)(P/p) c = (n + \)Kp\l", P c , (22) 
the second-order equation of hydrostatic equilibrium ([7]) , takes the dimensionless form of the Lane- Emden equation 

!(f 2 f) + ^»"=° ■ < 23 ' 

In terms of the dimensionless enthalpy #„(£) = H/H c , the dimensional included mass, mass density, average included 
mass density, and specific gravitational force are 

m(r) = ^ Pc o? ■ (~e0' n ), Pn(r) = Pc ■ 0£(O, Pn{r) := = Pc • (-3^/0, ff(r) := 47rp c a 2 (-^) (24) 

where prime designates the derivative ' := d/d£. The scale invariants are 

u:=-£0 n n /0' n , v n i=-&J9 ni u n :=(uv^ n -^=-e + ^e' n . (25) 

The Noether charge 

G n (o = ■ {e ■ Kf +^1)+^]} . ( 26 ) 

evolves radially according to 

rfG„ ,„ „. _ aH?\ /n-5\ iO /0i 2 (9" +1 



This non-conservation law expresses the radial evolution of energy density per mass shell, from entirely internal 
(6™ +1 /(n + 1)) at the center, to entirely gravitational (0^/2) at the stellar surface. 

Figure 1 shows the first integrals w n (u) for n — 0, 1, 2, 3, 4, 5. For n = 5, scaling is a variational symmetry so that 
(26) reduces to a conservation law for the Noether charge 

G 5 = ^-^ 2 K(^ + f ) + ^]=-^-(™ 5 3 ) 1/2 ' [-« 6 -u/3 + l]= constant . (28) 

For the Emden solution, W5 is finite at the stellar boundary u = 0, the constant vanishes, and 105(1*) = 5i>5 = |(3 — u) 
everywhere. 

For n < 5, v n diverges at the stellar radius £1 , but ui n — > oio n , a finite constant characterizing each Emden function. 
At the boundary u = 0, our density invariant w n {u) diverges as n[oU) 

(-eXJ^oWn-^ • (29) 

Table I lists these constants oUJ n , along with the global mass density ratios p c /p n (R) and the ensuing dimensional 
radius-mass relation M 1 ~ n = [(n + 1)K/G] n ■ (ow™" 1 /in)R 3 ~ n . Together with the well-known (HEi] third, fourth 
and fifth columns, all of this table follows directly from the regular solutions of the first-order equation (|21[) . In 
addition, the sixth and seventh columns express mass concentration in an original way. 



III. INCREASING POLYTROPIC INDEX AND MASS CONCENTRATION 

Emden functions are the normalized regular solutions of the Lane- Emden equation (|23p for which the mass density 
is finite at the origin, so that 6 n (Q) — 1, 6^(0) = 0. Each Emden function of index n is characterized by its first 
zero n (^i n ) — 0, at dimensionless boundary radius As an alternative measure of core concentration, we define 
the core radius £ corc implicitly by u(£, COTO ) :— 2, where gravitational and pressure gradient forces are maximal. This 
core radius, where w n 2 and the mass density has fallen to p nc0 rc/ Pnc ~ 0.4 for all polytropes n > 1, is marked 
by red dots in Figures 1, 2, 3. The sixth and seventh columns in Table I list dimensionless values for the fractional 



5 



TABLE I: Scaling Exponents, Core Parameters, Surface Parameters, and Mass-Radius Relations for Polytropes of Increasing 
Mass Concentration. Columns 3-5 are well-known [1, [H, |(| • Columns 6-7 present a new measure of core concentration. 



n 


&n 


fin 


pcn(R)/pn(R) 




Encore/ R — ^ncore/^1 


tic ore /A'if 


Radius-Mass Relation R 3 ~ n ~ M 1 ™/o^n 





-2 


2.449 


1 


0.333 


l 


l 


i? ~ M 1//3 ; mass uniformly distributed 


1 


dzOO 


3.142 


3.290 




0.66 


0.60 


R independent of 


1.5 


4 


3.654 


5.991 


132.4 


0.55 


0.51 


i? ~ M~ 1/3 


2 


2 


4.353 


11.403 


10.50 


0.41 


0.41 




3 


1 


6.897 


54.183 


2.018 


0.24 


0.31 


M independent of R 


4 


2/3 


14.972 


622.408 


0.729 


0.13 


0.24 




4.5 


4/7 


31.836 


6189.47 


0.394 


0.08 


0.22 




5 


1/2 


oo 


oo 








0.19 


R — oo for any M; mass infinitely concentrated 



core radius r ncoIe /R = £ nC ore/£i and fractional included mass m ncoie /M. Within the core u > 2, the internal energy 
dominates over the gravitational energy, so that for n > 1 , 

5 

w n (u) « io 5 (u) = -(3 -u) , 6»„(£) « 1 - £ 2 / 6 , for u n > 2, £ < £ core , (30) 

consistent with the universal density structure (8) all stars enjoy near their center. 
For n — 0, the mass is uniformly distributed, and the entire star is core. 

As < n < 5 increases, the radial distribution concentrates, and the envelope outside the core grows. With increas- 
ing core concentration: 

For 1 < n < 3, the radius R decreases with mass M. Nonrelativistic degenerate stars have n = 3/2. 

For n=3, the radius R is independent of mass M. This astrophysically important case is discussed in Section 
IV and the Appendix. 

For n > 3, the radius R increases with mass M. As n — > 5, the stellar radius increases £i„ — > 3(n + l)/(5 — n), 
the core radius shrinks £ core — » ^/10/3n, the fractional core radius r corc /R — £core/£i« —> 0.045(5 — n), 
m ncore /M -> 0.19, and w„ \/3j£vn -> 0. 



"O 

I 

II 




1.0 1.5 2.0 

u=dlog m/dlogr=3p/p 



FIG. 1: Dilution of polytrope mass density as the boundary is approached (it — > 0). All solutions are tangent to the same 
density structure w„(z) — > w$ — (5/3) (3 — u) at the center [u = 3), but differ for u < 2 outside the core. Approaching the outer 
boundary (u n — > 0), the density p n (r) falls rapidly, but uv„ := approaches a constant o^™ -1 so that w„ — > nlowJJ" 1 /!^] 1 '™ 

diverges, for n < 5. 
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TABLE II: Taylor Series and Picard Approximations #„pi c (£) to Emden Functions n (£) 



n 


Emden Function 6 n {£) and Taylor Series 


N n := 5/(3n - 5) 


Picard Approximation n Pic(£) := (1 + £ 2 /6N n )~ Nn 







-1 




1 


sin £/£ = 1 - £ 2 /6 + £ 4 /120 - C 6 /5040 + • • • 


-5/2 


(1 - £7i5) 5/2 = 1 - £ 2 /6 + £ 4 /120 - C 6 /10800 + • • ■ 


n 


1 - £ 2 /6 + n£ 4 /120 - n(8n - 5)/15120£ 6 + ■ • ■ 


5/(3n - 5) 


(1 + £ 2 /6N„)- N " = 1 - C 2 /6 + < 4 /120 - n(6n - 5)£ 6 /10800 + 


5 


(1 + C 2 /3)" 1/2 


1/2 


(l + £ 2 /3)- 1/2 



For n = 5, the mass is infinitely concentrated toward the center, and the stellar radius R = oo for any mass M. 
Scaling becomes a variational symmetry, so that the Noether charge G5 in (40) is constant with radius. 
For the regular solution this constant vanishes: 

G5 ~ [ ^T + ? ) + ^] = (w 5 3 ) 1/2 -("5-u/3-l) = , (31) 
so that = 1 — u/3, #5 — —^6*1/3. Integrating then yields 

5 (C) = (l+C 2 /3)- 1/2 , (32) 

after normalizing to 6*5(0) = 1. 
For n > 5 the central density diverges, so that the total mass M is infinite. 



IV. REGULAR EMDEN SOLUTIONS AND THEIR APPROXIMATIONS 



In place of u, we now introduce an equivalent homology invariant z := 3 — u = —d\ogp n /d\ogr, where p n := 
3m(r)/47rr 3 is the average mass density inside radius r. In term of z, w n , the characteristic differential equations (|20[) 
are 

dz d\ogw n d log to 

7^ \7 T = o T r = dlogr = ^ • ( 33 ) 

(6 — z){w n — z) 2 — z + w n /n 6 — z 

Incorporating the boundary condition, the first of equations (40) takes the form of a Volterra integral equation [9] 



w n {z) 



dz w n f Z + Wn/U \ « (5/J„)[l - (1 - z/3) J "] := w nPic (z) , J„ := (9n - 10)/(7 - n). (34) 
(3-z)(w„-z) 



The Picard approximation is defined by inserting the core values w n (z) rs (5/3)2 inside the preceding integral. For 
n = 0, 5, this Picard approximation is everywhere exact. For intermediate polytropic indices < n < 5, the Picard 
approximation breaks down approaching the boundary, where w n diverges as w n — > ii[quj™~ 1 /u] 1 ' n , and is poorest for 
n w 3. 

After obtaining w n (z) := —dlogp n /dlogr, either numerically or by Picard approximation, another integration 
gives [9[ 

e » - k-MMJ 1 * - {- jT „ K g:^L,) ) - p - =/ 3 » 5/2 " = *- < 3e » 

mW/i/ . (£)•/» ■ exp ^{ - sH " (i) 3 ^ (3?) 

-M/fi - f/fa. - ( . exp { jT d,{ (3 _ z)[ l M _ z] - 1}} . . (38) 

All the scale dependance now appears in the integration constants M and R(M), which except for n — 3 depends on 
M. Inserting the core values w n {z) (5/3)2 inside the integral, the Picard approximations 

e n Pic(0 = (l+f 2 /6^n)~ JV " , JV„:=5/(3n-5) (39) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

m/M 



FIG. 2: Normalized mass density profiles as a function of fractional included mass m/M, for polytropes of mass concentration 
increasing with n. The red dots mark the core radii, at which the densities stay near p(r COIC )/p c ~ 0.4, for all n > 1. For 
uniformly distributed mass (n = 0), the polytrope is all core. As the mass concentration increases (n — i- 5), the core shrinks to 
about 20% of the mass. 

to the Emden functions are obtained and tabulated in the last column of Table II. For polytropic indices n = 0, 5, 
this Picard form is exact. For intermediate polytropic indices < n < 5, the Picard approximation remains a 
good approximation through order £ 6 , but breaks down approaching the outer boundary. Unfortunately, the Picard 
approximation is poorest near n = 3, the astrophysically most important polytrope. 

Figure 4 compares three approximations to this most important Emden function, shown in yellow, whose Taylor 
series expansion is 

3 (£) = l - £ 2 /6 + £ 4 /40 _ (l9/5040)£ 6 + (619/1088640)£ 8 - (2743/39916800)£ 10 + ■ • • . (40) 
Tenth-order polynomial approximation to this Taylor series expansion 

1 - 0.1666667£ 2 + 0.025£ 4 - 0.0037698£ 6 + 0.0005686£ 8 - 0.00006872£ 10 , (41) 
shown in red, diverges badly for £ > 2.5 w 1.7 ^core- 
Picard approximation 

#3Pic(£) = (1 + 2£ 2 /15)- 5/4 = 1 - £ 2 /6 + £ 4 /40 - 13£ 6 /3600 + . . . ( 42 ) 
= 1 - 0.1666667£ 2 + 0.025£ 4 - 0.003611£ 6 + • • • , (43) 

shown in dashed green, converges and remains a good approximation over the bulk of the star, with < 10% 
error out to £ ~ 3.9, more than twice the core radius and more than half-way out to the stellar boundary at 
£13 = 6.897. This approximation suffices in white dwarf and ZAMS stars, except for their outer envelopes, 
which are never polytropic and contain little mass. Because it satisfies the central boundary condition, but not 
the outer boundary condition, the Picard approximation underestimates and overestimates outside 

£-3.9. 
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FIG. 3: Normalized mass density profiles as function of fractional radius r/R. The density is uniform for n — 0, but is 
maximally concentrated at finite radius for the n = 5 polytrope, which is unbounded (R = oo). The density at the core radius 
stays about p(r corc ) / p c ~ 0.4, for any n > 1. 



Pade rational approximation [Tol HT|: 

^3Pad = YJl^ywsl^yiQQg = 1-0.166667£ 2 +0.025£ 4 -0.00376984£ 6 +0.0005686£ 8 -0.0000857618£ 10 +- • • , 

(44) 

shown in dashed heavy black, is a simpler and much better approximation. By construction, it agrees with the 
series expansion (40) through fourth order. In fact, this Pade approximation is almost exact out to its first zero 
at £1 = 6.921, very close to the true outer boundary £ 13 = 6.897. 

These simple analytic approximations to #3(£) simplify structural modeling of massive white dwarfs and ZAMS stars. 



V. CONCLUSIONS 



We have explored how a symmetry of the equations of motion, but not of the action, reduces a second-order differ- 
ential equation to first-order, which can be integrated by quadrature. In scale-invariant hydrostatics, the symmetry 
of the equations yields a first integral, which is a first-order equation between scale invariants, and yields directly all 
the familiar properties of poly tropes. 

We observe that, like all stars, polytropes of index n share a common core density profile and defined a core radius 
outside of which their envelopes differ. The Emden functions n (^), solutions of the Lane-Emden equation that are 
regular at the origin, are finally obtained, along with useful approximations. 

The Appendix reviews the astrophysically most important n = 3 polytrope, describing relativistic white dwarf 
stars and zero age main sequence stars. While reviewing these well-known applications [a,|(|, we stress how these 
same mechanical structures differ thermodynamically and the usefulness of our original (Section IV) approximations 
to these Emden functions. 



Appendix: Astrophysical Applications of the n — 3 Polytrope 



The n — 3 polytrope, which is realized in white dwarfs of maximum mass and in the Eddington standard 
model for ZAMS stars just starting hydrogen burning, is distinguished by a unique M — R relation: the mass 
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FIG. 4: The exact Emden function #3 (f) (solid yellow) and its polynomial (red), Picard (green dashed) and Pade (heavy black 
dashed) approximations. Even in this worst case, the Picard approximation holds out to twice the core radius at 2^3 Corc = 3.3, 
before breaking down near the boundary. The Pade approximation is indistinguishable from the exact solution, vanishing at 
£1 = 6.921, very close to the true boundary at £13 = 6.897. 



M = Att{qlo^) (K/ttG) 3 / 2 is independent of radius R, but depends on the constant K := P/p 4 / 3 . In these stars, the 
gravitational and internal energies cancel, making the total energy W — fl + U = 0. Because these stars are in neutral 
mechanical equilibrium at any radius, they can expand or contract homologously. 



1. Relativistic Degenerate Stars: K Fixed by Fundamental Constants 



The most massive white dwarfs are supported by the degeneracy pressure of relativistic electrons, with number 
density n e = p/p, e mH, where ran is the atomic mass unit and the number of electrons per atom p e = Z/A = 2, 
because these white dwarfs are composed of pure He or C 12 /0 16 mixtures. Thus, Ky/D = (/ic/8)[3/7r] 1 / 3 TOij/i e -4 / 3 
depends only on fundamental constants. This universal value of KyyD leads to the limiting Chandrasekhar mass 
A/ C h = (ir 2 /8Vl5)M*///2 = 5.824M Q / M 2 = 1.456M Q • (2//i) 2 



2. Zero-Age Main Sequence Stars: Mass and K(M) Dependent on Specific Radiation Entropy 



In an ideal gas supported by both gas pressure P gas = IZpT / ji := f3P and radiation pressure P ra d = aT 4 /3 := 



(1 — /3) P, the radiation/gas pressure ratio is 

Pad >= = r 3 W 

Pgas ' P P ' in 

The specific radiation and ideal monatomic gas entropies are 



i' ad 



4aT 3 



S ga8 (r)=(^)-log 



so that the gas entropy gradient 



rT(r) 5 / 2 



p(r) 



dlogP \2p, 



(A.45) 



(A.46) 



(A.47) 
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depends on the difference between the adiabatic gradient V a d = 2/5 and the star's actual thermal gradient V 
dlogT/dlogP, which depends on the radiation transport. 

Bound in a polytrope of order n , the ideal gas thermal gradient and gas entropy gradient are 



V:=l/(n + l) , -J^= - • — — --1 . (A.48) 
dlogP V /i / L2(n + 1) J 

For n > 3/2, the thermal gradient is subadiabatic, the star's entropy increases outwards, so that the star is stable 
against convection. 

ZAMS stars, with mass OAMq < M < 15OA/0, have nearly constant radiation entropy S T3 d(M), because radiation 

transport leaves the luminosity generated by interior nuclear burning everywhere proportional to the local transparency 
(inverse opacity) k" 1 . Assuming constant S ra ^(M), we have Eddington's standard model, an n — 3 polytrope with 
S ra d(M) = 4(ft//i) • (1 - and 

K{M) = P/pW = {[3(1 - /3)/a]CW)Y /3 , (A.49) 
depends only on (3(M), which is itself determined by Eddington's quartic equation [1, [H, |(| 

1-/3 /M/i 2 \2 3^10 w 3 / he \3/2 
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The luminosity 

L = L Edd [l - p(M)} = L Edd ■ (O.OO3) M 4 /?(M) 4 (M/M ) 3 , (A.51) 

depends on the Eddington luminosity L Edd :— AttcGM/ k p through the photospheric opacity K p . This mass-luminosity 
relation is confirmed in ZAMS stars: on the lower-mass ZAMS, /3 w 1, L ~ A/ 3 ; on the upper-mass ZAMS, /3 ~ 
(M^ 2 /M*)- 2 <1, L~M|, 
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